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1. Introduction 

Finding or excluding the Higgs boson of the Standard Model (SM) of elementary particle 
interactions is amongst the top priorities of the CERN Large Hadron Collider (LHC) 
experiment. So far the Higgs has escaped direct detection, although the year 2011 run 
culminated in intriguing hints pointing towards a Higgs boson in the ^ 125 GeV mass 
range [1, 2, 3, 4]. 

In simple extensions of SM, the Higgs can be easily hidden from the standard searches, 
explaining the absence of a clear Higgs signal in early LHC data. In a typical scenario, 
the Higgs boson decays dominantly into a pair of weakly interacting stable particles, which 
then escape the detector and are only seen as missing energy. If the production of the 
Higgs boson is unaffected by these new physics degrees of freedom, the invisible decays 
effectively reduce the cross section of the final states that are looked for in the standard 
Higgs searches, since only a subleading part of the total number of Higgs bosons produced 
end up in these final states. In this case the exclusion limits quoted by the experiments do 
not apply as such, and a light Higgs boson could still be there. By adjusting the relative 
decay width of the invisible channel, some of the excess events in Higgs boson searches can 
also be explained [5] . 

The coupling of the Higgs boson contributes to the mass of the particle it couples to, 
and hence this new weakly interacting state must be massive in order to dominate the other 



decay channels. By construction it is then a weakly interacting massive particle (WIMP), 
and could perhaps also provide plausible dark matter candidate. This scenario has been 
considered recently in the case of a scalar singlet dark matter [6] and vector dark matter 
[7]. The case of a singlet fermion was studied in [8] and [9]. 

In this paper we point out that this construction may arise as a consequence of ad- 
dressing the naturality in the Higgs sector. A concrete example is provided by the minimal 
walking technicolor (MWT) model [10] where a fourth generation of leptons is required to 
cancel the global anomaly associated with the strongly interacting sector responsible for 
electroweak symmetry breaking. 

More generally, we consider an effective theory for a heavy fourth generation lepton 
doublet couphng to an effective SM hke (possibly composite) Higgs field. We investigate 
if it is possible for the neutrino to provide the dominant decay channel of the (in this case 
composite) Higgs boson, and simultaneously produce the right amount of relic density to 
contribute the observed dark matter abundance. This scenario is significantly different 
from the SM with four sequential generations (SM4) model, where a complete generation 
of quarks and leptons is added to the SM. In the SM4 model the loop induced Higgs 
coupling to gluons is enhanced by the two new heavy quarks running in the loop. This 
effect significantly enhances the production cross section of the Higgs boson in the gluon 
fusion process, and effectively negates the effect that the fourth generation neutrino could 
have in hiding the Higgs boson. In the MWT model there are no fourth generation quarks, 
so the production cross section of the Higgs is unchanged, but the decays may be strongly 
affected by the new neutrino. For previous results on Higgs decaying into a pair of neutrinos 
see [11]. Our main finding is that achieving both the observed dark matter relic density 
and strong enough invisible decay width is heavily constrained in this type of model. 

In section 2, we will first review the aspects of the MWT model which are relevant for 
the scenario we have outlined in this section. Then, in section 3 we present our numerical 
results and constraints in light of present data from colliders and dark matter experiments. 
In section 4 we present our conclusions and outlook. 

2. Minimal Walking Technicolor and the Fourth Generation of Leptons 

2.1 Model Lagrangian and mass terms 

Since we imagine the Higgs to result from a gauge theory confining at the electroweak 
scale, the low energy degrees of freedom are described better in terms of a chiral effective 
theory than using the fundamental techniquark and -gluon degrees of freedom. The form 
of the effective theory is fixed by the underlying chiral symmetry breaking pattern. The 
simplest possibility is SU(2)LxSU(2)i^ ^SU(2)y, which has three Goldstone bosons. By 
coupling with the electroweak currents these are absorbed into the longitudinal degrees of 
freedom of the weak gauge bosons by the Higgs mechanism, and only the CP-even scalar 
Higgs remains in the physical spectrum to couple with the matter fermions. 

We denote the left-handed fourth generation lepton doublet by Ll — (N^, El) and the 
right-handed SUl(2) singlets as Er and Nr. The interactions between the scalar sector 



and these leptons up to and including dimension five operators are given by 

-CMass = {vLlHEr + h.C.) + CdLlHNr 

+ ^{LYH){H^Ll) + ^{H^HWrNr + h.C. (2.1) 

where H — ir^W and A is a suppression factor related to the more complete ultraviolet 
theory like extended Technicolor (ETC) providing a more microscopic origin for these 
interactions. The first terms in Eq. (2.1) lead to the usual (Dirac) mass for the charged 
fourth generation lepton, and the remaining terms allow for more general mass structure 
of the fourth neutrino. After symmetry breaking the effective Lagrangian (2.1) gives rise 
to a neutrino mass term: 

1 _^ ( Ml mD\ ^ , .^^. 

2 yrriD Mr J 

where ul = {Nl^N^)^ ^ mo = Cd'^/V^ and Ml^r = C^^rv^ j^k^ where v is the vacuum 
expectation value of the effective Higgs field. The special cases are a pure Dirac and a pure 
Majorana neutrino which are obtained, respectively, by discarding dimension five operators 
and by removing the right handed field Nr. The most general mass matrix contains, even 
after the field redefinitions, one complex phase. However, in this paper we shall restrict 
ourselves to the case of real mass matrix. Mass eigenstates are two Majorana neutrinos 
which are related to the gauge eigenstates by a transformation 

AT = OriL + pO^nl , (2.3) 

where A^ = (A^i, N^)^ and O is an orthogonal 2x2 rotation matrix, where the associated 
mixing angle is 

tan 2^=—^^^. (2.4) 

Mr -Ml ^ ' 

The phase-rotation matrix p = diag(/?i,/92) is included above to ensure that the physical 
masses mi_2 are positive definite. Indeed, the eigenvalues of the mass matrix in (2.2) are 

A± = ^(^Ml + Mr± ^{Ml - MrY + Amiy (2.5) 

Because the signs and relative magnitudes of Ml^r and ttid are arbitrary, the eigenvalues 
\± can be either positive or negative. However, choosing independent phases as p± — 
sgn(A±) we get positive m± — \\±\ as required. For our purposes it will be convenient 
to express everything in terms of the physical mass eigenvalues rui > m2 and the mixing 
angle sin^ instead of the Lagrangian parameters M^,, Mr and tud- While working with 
physical parameters has obvious advantages, the downside is that the connection between 
the physical and the Lagrangian parameters is not always straightforward. 

The role of the phase-rotation parameters pi and p2 has been discussed in detail in 
[12]. Here we simply point out the feature that the physically relevant parameter is the 
relative sign of pi and p2^ 

Pip2 = Pi2 = ±1, (2.6) 



which divides the parameter space into two parts. The typical hmits of purely left- or right- 
handed neutrinos and the Dirac-limit are contained in the pi2 — —1 part of the parameter 
space. We will give our results for both values of this parameter. 

Naturally also the mixing angle can have positive or negative values depending on the 
Lagrangian parameters M^,, Mr and ttid- However, our results depend only weakly on the 
sign of the mixing angle and thus we will present our results only for the positive mixing 
angle values. 

Although the relevant degrees of freedom and their interactions described above are 
generic and may arise in different beyond the Standard Model scenarios, it is good to have 
at least one particular microscopic realization available. We consider the MWT model, 
where the electroweak symmetry breaking is driven by the gauge dynamics of two Dirac 
fermions in the adjoint representation of SUtc(2) gauge theory. The key feature of this 
model is, that it is (quasi) conformal with just one doublet of technifermions [10]. This 
feature is essential for a technicolor theory not to be at odds with electroweak precision 
measurements. However, since the technicolor representation is three dimensional, the 
number of weak doublets is odd and hence anomalous [13]. A simple way to cure this 
anomaly is to introduce one new weak doublet, singlet under technicolor and QCD color 
[10, 14] in order not to spoil the walking behavior and to keep the contributions to the 
oblique corrections as small as possible. Hence, the model requires the existence of a fourth 
generation of leptons. The anomaly free hypercharge assignments for the new degrees of 
freedom have been presented in detail in [14]; for the techniquarks we have 

y(Qi) = |, Y{UR,DR)^(^-±^y-^^, (2.7) 

while for the fourth generation leptons we have 

In the above equations y is any real number. Choosing y — 1/3 makes the techniquarks 
and the new lepton doublet appear exactly as a regular standard model family from the 
weak interactions point of view. The heavy fourth generation neutrino becomes a natural 
dark matter candidate provided that it is stable; to ensure the stability one can postulate 
a discrete symmetry. We note that under this hypercharge assignment fractionally charged 
techniquark-techniquark or techniquark-technigluon bound states may form in the early 
universe. The existence of such charged relics is severely constrained ^ . However, according 
to Refs. [15] and [16] there seems to exist an open window for relics with masses ^ 0.1- 
10 TeV. This is exactly the range where one would expect the mass of the technihadron 
states to lay, as the natural scale is of (9(TeV). Furthermore, when interpreting especially 
the collider limits, one should keep in mind that, as the relic here is a composite state 
its production and signals in colliders differ from those of charged elementary particles. 
Additional limits for these relics can follow from WMAP data as discussed in [17]. However, 



^Related to this specific bound state of tecfiniquark and tecfinigluon tfiere is a recent discussion in [21]. 



Ref. [17] do not plot their exclusion region to the charge values of our interest. One should 
also notice, that the WMAP data is always analyzed using specific cosmological model, thus 
these constraints are model dependent and should not be taken too strictly. Moreover, we 
do not know for certain if these bounds states would be even created in the first place 
in the early universe. To give specific constraints for their fate, more complete analysis 
should be performed e.g. including determination of the interactions of this state with SM 
particles and the calculation of the relic density. Because this state in the end is part of 
the other sector (TC) of the underlying model, we postpone this study to future work, and 
concentrate here to the case where we assume that either this bound state does not exist 
or that it does not affect the cosmology^. 

On the other hand one should notice, that if the hypercharge assignment is changed to 
^ = 1 there are no fractionally charged states, and the fourth generation leptons are doubly 
and singly charged and their contribution to the dark matter relic density is disfavored. 
However, the techniquark-technigluon bound state containing L)-techniquark becomes elec- 
trically neutral. This state will then be the DM candidate instead of our new neutrino. 
This kind of model has been studied previously in [22]. In this case our analysis and results 
presented here can be directly interpret as an analysis for the techniquark-technigluon dark 
matter. A somewhat similar scenario of a composite DM particle from a strongly interact- 
ing hidden sector has been considered in [23]. As yet another model building alternative, 
instead of single lepton doublet with hypercharge Y(Ll) = —3/2, one can saturate the 
Witten anomaly by introducing three SM-like lepton doublets, and a stable heavy neutrino 
among these three becomes a plausible dark matter candidate. Various further possibilities 
of nonsequential generations beyond the SM have been considered in literature [24, 25]. 

Finally we note that in MWT, the global symmetry breaking pattern is SU(4) -^ 
SO (4), with nine Goldstone bosons. Three of these are absorbed into the longitudinal 
degrees of freedom of the weak gauge bosons, and the low energy spectrum is expected 
to contain six quasi Goldstone bosons which receive mass through extended technicolor 
interactions [26, 27, 28]. Their phenomenology has been investigated elsewhere [29, 30, 21]. 
For the hypercharge assignements we consider, the conservation of hypercharge allows 
only the effective SM-like scalar Higgs to couple to the fermions. Hence, to consider the 
interactions between fermions and the scalar sector, the interactions introduced in Eq. 
(2.1) are sufficient. 

Of course it is possible that the (composite) Higgs is not the only source of the fermion 
masses, but there are other (composite or even fundamental) scalars whose condensation 
leads to mass terms for the matter fields. To illustrate such possibilities, we consider as an 
alternative to the model Lagrangian (2.1), the case where the right-handed neutrino mass 
originates from a Standard Model singlet scalar field S. 



^Ma«s = (vLlHEr + h.C.) + CdLlHNr 
Ci 

A 



+ % {Ll'H) (H^Ll) + CrSN^^Nr + h.C. (2.9) 



^For positive effects of charged relics in cosmology see e.g. [18] 



This model is similar to the usual see-saw neutrino mass generation mechanism, although 
here the singlet S does not need to be a fundamental scalar. To specify the model com- 
pletely one should give a potential for S. However, none of the parameters of this potential 
are needed in our analysis since we may assume that the vacuum expectation value for S 
is generated through interactions with the Higgs, i.e. we do not need additional sources of 
spontaneous symmetry breaking for the dynamics of the S-field. In what follows, we will 
refer to the scenario with just only the doublet Higgs field as Scenario I and to the case 
with Higgs and a singlet scalar as Scenario H. 

2.2 Couplings 

For the analysis of the Higgs decay branching ratios and relic density we need the cou- 
plings of the neutrino mass eigenstates to the Higgs boson and to the weak gauge bosons. 
These are easily found out by applying the appropriate phase- and rotation transforma- 
tions defined in the previous section. We shall write down only the terms relevant for our 
calculations. For the Z and W^ bosons we find that 

W+Nl-z^'El = sine W+N2L-f^EL + ■■■ 
Z^Nli'^Nl = sin2 Z^N2Ll''N2L 

+ \ sin 29 Z^ {NiLl^N2L + N2L1''Nil) + ■■■ , (2.10) 

where the omitted terms contain interactions of the heavy Ni field only. These couplings 
are diagonal in the mixing and therefore do not involve the phase factor pi2. However, 
neutral current involves mixing and these couplings do depend on pi2. One finds: 

N2-f^Z^PLNi + Ni-f^Z^PLN2 = MP + 075)7^ ^M^i, (2.11) 

where 



a 



U1 + P12) and /3=|(l-pi2). (2.12) 



Thus, for P12 = — 1 the neutral current interaction of our WIMP is purely axial vector 
and for pi2 = +1 purely vector. Usually in the literature dealing with the interactions of 
Majorana neutrinos, only the first possibility is mentioned, although e.g. [9] considers both 
operators. 

The effective interaction terms involving the Higgs and the lighter neutrino eigenstate 
are 

Cnh = ^ ( C^2hN2N2 + C^^hNiia - /3j^)N2 

+ ^-C^^h'N2N2) + ^h' + ---, (2.13) 

where we have again omitted the interaction terms which do not contain N2 and hence are 
not needed in our analysis. 

The interactions between the Higgs and the neutrino can be generically described by 
the Lagrangian (2.13) for scenarios I and II we have introduced. The factors a and /3 are 
defined in Eq. (2.12) and the factors C22, C21 and C22 are given in Table 1. 





Scenario I 


Scenario II 


^22 


1- \^ii?2eR- 
-\pi2^\nA9 R- 


sin^ 9 

lpi2sm29R^ 

lsm^9{l-co8^9R-) 



Table 1: Coefficients of the Lagrangian (2.13) for two the 
distinct mass generating scenarios described by Eqs. (2.1) 
and (2.9). We have defined R± = l±pi2^. 



2.3 Oblique constraints 

The fourth generation of leptons is 
constrained by current accelerator 
data. From LEP we know that the 
charged lepton E has to be more 
massive than the Z boson and if 
the fourth generation neutrino has 
standard model interaction strength, 
it needs to be heavier than Mz/2 
in order to evade the constraint from Z-pole observables. In the case of neutrino mixing 
considered in this work, the lighter state can have a substantial right-handed component 
and hence interact only very weakly. This could allow this state to escape the LEP bounds 
even when its mass is less than Mz/2 (see e.g. [12]), but in this work we will limit to 
the case 1712 > Mz/2. In addition to these direct bounds, the parameters of the fourth 
generation leptons are constrained by oblique corrections, i.e. due to their contribution 
to the vacuum polarizations of the electroweak gauge bosons. These contributions are 
conveniently represented by the S and T parameters [31]. 

The oblique corrections in MWT model with the general mass and mixing patterns 
considered here have been studied in detail in [32] and [55]. We also note that there 
exists two extensive fits performed by the LEP Electroweak Working Group (LEPEWWG) 
[33] and independently by the PDG [34]. Both fits find that the SM, defined to lie at 
(5, T) = (0, 0) with mt = 170.9 GeV and tuh = 117 GeV, is within la of the central value 
of the fit. The two fits disagree slightly on the central best-fit value: LEPEWWG finds 
a central value (S, T) = (0.04, 0.08) while including the low energy data the PDG finds 
(5, T) = (0.03,0.07). Since the actual level of coincidence inferred from these fits depends 
on the precise nature of the fit, we allow a broader range of S and T values, roughly 
corresponding to the 3a contour. From the results of [32] it can be inferred that these 
values can be accommodated easily within the parameter space of the leptonic sector. 

For reference we show the allowed mass spectrum of the fourth generation in figure 1 
for a choise of the lightest neutrino mass 1712 = 62 GeV, and sin^ = 0.1,0.2 or 0.3. As 
was discussed in [32], the limiting factor is the T-parameter, which fixes the ratio of the 
two masses mi and ttie to a narrow range. This ratio is slightly dependent on the value of 
the mixing angle, as can be seen in figure 1. Here we have used pu = — 1. For pi2 = +1 
the T-parameter is generally a bit smaller and thus the allowed range of masses is slightly 
wider. The absolute scale of the masses affects the S parameter slightly, but the constraint 
from S is much weaker than the one from T. 

Here we are interested in finding parts of the parameter space, where the precision 
constraints are met, the lightest neutrino mass eigenstate provides the correct relic density 
to match the observed DM abundance, and the invisible decay channel H -^ N2N2 is 
the dominant decay channel of the composite Higgs boson. Our strategy for scanning the 
parameter space is as follows: We scan the two-dimensional parameter space defined by 
the mass of the lightest neutrino mass eigenstate 7712 and the neutrino mixing angle sin0. 
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Figure 1: The allowed mass spectrum of the fourth generation leptons for m2 = 62 GeV, p\2 = — 1 
andsinl9 = 0.1,0,2 or 0.3. 



For each point in this plane we scan over the subspace defined by the mass of the heavier 
neutrino vn\ and the mass of the charged lepton tti^, and select the point in this plane that 
gives the most suitable value for the oblique parameters S and T. In the case p\2 = +1, 
practically the whole (7712, sin0)-plane is allowed in terms of the oblique constraints. That 
is, for every point in that plane there is a configuration of the values of vn\ and vfiE that 
produces acceptable values for S and T. If p\2 = — 1 the T parameter gets larger with 
large values of sin^ and the values above sin^ > 0.45 are ruled out. 

We then calculate the invisible decay width of the composite Higgs boson and the relic 
density of the light neutrino in this point of the parameter space. In the following sections 
we will describe the evaluation of the relic density, the constraints from earth-based direct 
dark matter searches, and the calculation of the invisible decay width of the Higgs boson. 
We will then present the results in section 3. 

2.4 Higgs Decay 

The Higgs coupling to the lightest neutrino is given in equation (2.13). The resulting 
tree-level decay width is 



GFim^C^^fma / /2m2 

I H^N2N2 = -^ \ ^ ~ [ 




The invisible branching ratio of the composite Higgs boson is defined as 

^H^N2N2 



Rr 



SM' 



^H^N2N2 + ^H 



(2.14) 



(2.15) 



where F^ is the total Higgs decay width in the SM. The relative decay width to a given 
SM decay channel H -^ XX is then modified by a factor of (1 — Rv)'- 

Rxx = p ^"""-TrSM = (1 - Rr)RTx, (2-16) 

where R^^ = ^H^xx/^^ is the corresponding branching ratio in the SM. Since the 
production cross section of the Higgs boson in our model is equal to the SM, this suppression 
effectively results in a suppression of the total cross section for a given Higgs boson search 
channel. 

2.5 Relic Density 

Here we will summarize and update the relic density analysis which is originally given for 
this particular model in work [12]. We calculate the DM abundance ^n^ in the standard 
way, using the Lee- Weinberg equation [35] for scaled WIMP number density: 

^=<^=M^(/(.)-/JW), (2.17) 

where we have introduced the variables 

/(x) = ^^, and x = ^^, (2.18) 

Se 1712 

where m2 is the WIMP mass and se{T) is the thermal entropy density at the temperature 
T. H{T) = (87rp(T)/3M2j)V2 is the Hubble parameter and (av) is the average WIMP 
annihilation rate which expression is defined below in Eq. (2.20). We assume the standard 
adiabatic expansion law for the universe and use the standard thermal integral expressions 
for Se and for H{T). Freeze-out temperature for our WIMPs is typically T ^ (9(1 — 
10) GeV and thus the uncertainties in se related to the QCD phase transition do not 
affect our analysis. The present ratio of A^2-number-density to the entropy density /(O) is 
solved numerically from Eq. (2.17) after which the fractional density parameter 0.^2 of the 
Major ana WIMPs follows from 

f^jv, c:^ 5.5x10" ;g|/(0). (2.19) 

From Eq. (2.17) one sees that the rehc density /(O) essentially depends on the ratio {av)/H; 
the smaller the ratio, the less time the WIMPs can remain in thermal equilibrium and thus 
the larger is their relic density. It can be shown that the dependence is in fact almost 
linear: ^^2 ^ H/ {av) (see e.g. [36]). As we assume the standard expansion history of the 
universe so that the H is known, the solution /(O) is determined by the annihilation cross 
section {crv). 

For the thermally averaged annihilation cross section we use the expression [37]: 

1 f^ fs 

<""> = 8.n|rAl(=y) /„, '"^'^' - ""'D^Vf )..o.{.) (2^20) 



where Ki{y)s are modified Bessel functions of the second kind and s is the Mandelstam 
invariant. For the total cross section atot we considered the N2N2 annihilation to the final 
states including all open fermion, gauge boson and scalar channels 

N2N2 -^ ff, W^W-, ZZ, ZH^ and H^H^ , (2.21) 

Here H^ is the effective light "SM-like" Higgs state appearing in the mass operators (2.1) 
and (2.9). We did not take into account the WIMP annihilations to technifermions. This 
is because the technifermions would give only a small contribution to the fermionic annihi- 
lation channel in particularly in the case of heavier WIMPs, of which we are not interested 
here. Also in the WIMP mass Scenario II, the annihilations to scalars S were omitted 
assuming that these scalars are heavy. The cross sections for each channel shown in (2.21) 
were calculated without further approximations and all 5-integrals were solved numerically. 
For these computations the needed WIMP-Higgs and WIMP-gauge bosons couplings were 
given in section 2.2. We also assumed that the unstable heavier neutrino state A^i has 
decayed before the N2 freeze-out. Thus, the particles present during the freeze out are just 
the Standard Model particles and the annihilating WIMP. 

Let us mention that the WIMP annihilation cross section to fermionic final states used 
in this work differs slightly from the one used in [12] (Eq. 3.6). This is because here we 
took into account the full expression of the Z-boson propagator (in unitary gauge) when 
calculating this cross section, as in [12] only the part of the propagator proportional to ^^^y 
was used. However, the new terms, following from the use of the full Z propagator, have 
only minor impact on the cross section and thus on the final WIMP density in the WIMP 
mass range of our interest. For other applications the new terms can be relevant and thus 
we give the corrected cross section in the Appendix. More details about the computation 
of the cross sections can be found from [12]. 

Before going to the results let us summarize shortly what parameters affect on the 
relic density: the relic density is controlled by the annihilation cross section, which scale is 
set by the mixing angle sin 6 and the mass of the WIMP 777,2. The Higgs-boson mass tuh 
and the WIMP mass scenario affects considerably to the relic density for WIMP masses 
777,2 ^ 777,if/2. Also the phase pu and the mass 777,1 have an impact, mostly through the 
Higgs couplings, on our results. However, the mass of the charged lepton 777,^ affects only 
very weakly on the relic density analysis. A more detailed characterization of the effects 
of different parameters on our results is given in the next section. 

3. Results 

3.1 Invisible Decay Width 

As discussed in section 2.3 we do a scan over the (777,2, sin 0)-plane and choose suitable 
values for 777,1 and 777,^ in each point of the plane, to meet the (5, T)-constraints. We then 
calculate Rr and the relic density for each point. The decay width of the Higgs is com- 
puted at tree level, except for the gg-channels where the leading logarithmic corrections 
are taken into account. The results are shown in figures 2, 3 and 4 for Higgs masses of 120 
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Figure 2: The invisible branching ratio Rr of the composite Higgs boson for uih = 120 GeV in 
scenario I (upper panel) and scenario II (lower panel). pi2 = — 1(+1) in the left (right) panel. The 
black dots show points of the parameter space that produce acceptable relic density. The area above 
the white dashed line is ruled out by the XENONlO-results, and the continuous white line shows 
the XENONlOO-limit, as discussed in section 3.2. We do not show points of acceptable relic density 
above sin^ > 0.42, since these are in any case ruled out by the direct detection experiments. In the 
left panels the area above the topmost solid line (yellow), is ruled out by the oblique constraints. 



GeV, 130 GeV and 145 GeV, respectively. The colormap shows the value of the invisible 
branching ratio Rr, and the black dots show the points where the neutrino relic density 
has the correct value. The area above the solid white line is ruled out by the earth based 
direct detection dark matter searches. 

We shall now characterize the relic density results shown in Figs. 2-4 more carefully. 
As the WIMP mass 7712 and the mixing angle sin 9 set the scale for (av) these parameters 
also have the strongest impact on the final WIMP abundance ri^Vs- For this reason we 
have projected all the suitable model parameter sets, which produce the measured DM 
density Qdm - 0.19 - 0.23 [34] (consistent with combined WMAP7+H0 results [38]), in 
the (7712, sin 0) plane. Our results are also sensitive on the mass of the light composite 
Higgs particle ttih, especially in the WIMP mass region 7712 ^ tuh/'^, where the WIMP 
annihilation cross section gets enhanced due to the increase of Higgs s-channel process at 
the Higgs pole. Our results also depend on the phase pi2 which can be seen from the 
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Figure 3: Same as figure 2 but for ttih = 130 GeV. 

WIMP-Higgs interaction terms in Table 1. Finally, the mass of the new charged lepton 
state, rriE, has only a subleading effect on the annihilation cross section and hence on 
the relic density. Indeed, the charged state contributes to the cross section only as a 
virtual state in the t-channel process in WIMP annihilations into VF+VF" final states. In 
particular, for WIMP masses ni2 < Mw of which we are interested here, this annihilation 
channel has only very small effect in the cross section integral {av) given in Eq. (2.20). In 
principle similar arguments, for the subleading effects in annihilation cross section, holds 
also for the heavier neutral state A^i with mass rui, as also A^i contributes only as virtual 
state in the t-channel annihilations into ZH^ and H^H^ final states. However, for WIMP 
mass Scenario I, as the Higgs couplings given in Table 1 include the mass ratio of the 
neutral states 7711/7712, also 7711 will affect the annihilation cross section and thus the final 
WIMP abundance for the WIMP masses 7712 ^ 77iif /2. 

For WIMP masses 7712 < Mw the WIMP annihilation cross section is determined by 
the WIMP annihilations to SM fermions via Z-boson in 5-channel process. The annihilation 
cross section in this case is directly proportional to mixing angle: (av) oc sin^ 9 (this can 
be immediately realized from the Z-WIMP interaction term given in Eq. (2.10)). Now for 
a fixed mixing angle the cross sections gets enhanced at the Z-pole. As the DM density 
is inversely proportional to the annihilation cross section, the DM density will in principle 
decreases in this case. Thus, to keep the ^n^ constant, indicating that the cross section is 
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Figure 4: Same as figure 2 but for ttih = 145 GeV. 

constant, we need to decrease the value of the mixing angle to compensate for the effect 
of the Z-pole. This produces a dip on the suitable mixing angle values for WIMP masses 
7712 ^ Mz/2. 

Similarly the effect of the Higgs resonance can be seen in Figs. 2 -4 as a dip in the 
mixing angle values for the WIMP masses 7712 -^ 777/^/2. For the WIMP mass Scenario 
I the effect of the Higgs pole is more dramatic than the effect of the Z-pole. This is 
so because in the Higgs case the cross section depends only weakly on the mixing angle 
and thus the mixing angle is not able to sufficiently compensate for the pole. Indeed, for 
WIMP mass values 7772 ^ tuh/'^ the mixing angle suppression is insufficient to keep the 
cross section small enough for production of the correct relic density. Thus the dip in the 
suitable mixing angle values for the WIMP masses 7772 -^ ^^if /2 becomes very deep and 
eventually it cuts the suitable parameter space well before 7772 reaches the value 777/^/2. 
However, for Scenario II the cross section is again proportional to the mixing angle, as can 
be realized from the WIMP-Higgs coupling factor C22 given in Table 1, and thus the Higgs 
pole can be compensated for by decreasing the mixing angle and the correct relic density 
can be obtained for larger range of 7772, sin 9 parameter points. Let us also mention that in 
Scenario II the similarity of the DM density curves in different pi2 cases follows from the 
fact that, as the C22 factor is independent of the pu parameter in this case, the (av) for 
light WIMP masses (7772 < 777^) is effectively independent of the pu factor. 
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In all our plots we have set an upper limit for the mixing angle sinO < 0.42. This 
cuts the suitable DM density regions/curves for heavy WIMP masses (7712 ^ ttih/'^?) e.g. 
in lower panels of Fig. 2 DM density curve is cut for WIMP masses larger than ^ 60 
GeV. Similar cuts can be realized in other figures. There are two reasons for setting this 
cut/upper bound. Firstly, we are interested on relatively light WIMPs i.e. 45 GeV < 7712 < 
75 GeV, and in this mass region sin0 values larger than ^ 0.4 are excluded already by the 
null results of direct DM search XENONIO experiment, as is explained in Sec. 3.2 and 
shown in the Figs. 2 -4 as the white dashed line. Thus, we do not need to plot the curves 
for heavier DM masses. Moreover, the model results for heavier WIMP masses were already 
shown in [12]. Secondly, in this work we are particularly interested in the Higgs decays into 
invisible DM channel, and since this decay is kinematically allowed only for 7712 ^ TTi/f /2, 
in this case plotting the DM density curve for much larger DM mass values than ^ m/f /2 
is not interesting. Saying all this, one should still keep in mind that in principle the DM 
density curve continues in all figures to the right (as shown in [12]), and that for large 
DM masses there still might be suitable parameters space left i.e. (7712, sin 0) points which 
produce correct DM density and are consistent with the EW precision measurements and 
are not excluded by the DM search experiments. 

3.2 Cryogenic limits 

Here we summarize and update the strongest observational constraints for our WIMP 
candidate following from the direct DM detection experiments XENONlO-100 [39, 40, 48]. 
Part of these constraints were already shown in the previous result section, and here we 
demonstrate how these constraints were derived and also analyze the constraints set by 
these experiments in more detail. Other weaker constraints arising from LEP measurement 
of Z-decay width and from CDMS [41] were studied in [12], and constraint coming from 
indirect DM detection in particularly from neutrino detectors like Super-Kamiokande were 
studied in detail [42]; see also indicative limits given in [12] ^. Further, the limits following 
from FERMI-LAT gamma ray data for this model will be studied elsewhere [47] . 

We can set constraints for our model using both spin- dependent (SD) and spin-independent 
(SI) WIMP-nucleon interactions limits given by XENONlO-0 collaborations. For our model 
the spin-dependent WIMP-nucleon interactions proceed via Z-boson exchange, while the 
spin-independent interaction is Higgs mediated. As demonstrated in [12], for Higgs masses 
rriH > 200 GeV, constraints arising from spin-dependent WIMP-nucleon cross section lim- 
its were more stringent than those following from the spin-independent cross section limits. 



^IceCube and SuperKamiokande collaborations have updated their constraints for spin-dependent 
WIMP-nucleon cross sections in [45] and [46] respectively. When comparing these new constraints with 
the older SUSY WIMP constraints given in [43] and [44] by Super-Kamionkande and IceCube collabora- 
tions, respectively, one notices that the WIMP-nucleon cross section especially for heavier WIMP masses 
(m ^100 GeV) is now more constrained. These constraints are model dependent and are given for SUSY 
models. Thus they can be considered only as indicative limits for our model. In the analysis done for our 
model in [42] and [12], the same Super-Kamiokande data as in the older SUSY analysis in [43] has been 
used. As an outcome, only WIMP masses between ^ 100 — 200 GeV in our case were excluded, and one 
can conclude that in the WIMP mass range {^ 45 — 75) GeV, which we consider here, these new limits are 
less severe for our model when compared to the XENONIOO limits given below. 
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However, here we consider lighter Higgs masses and as we are at the same time interested on 
relatively light WIMPs, the constraints arising from the spin- independent WIMP-nucleon 
interactions become stronger than the ones following from spin-dependent WIMP-nucleon 
interactions. It actually turns out that the WIMP mass Scenario I is basically excluded 
by the spin- independent limits (in the case of light Higgs masses). This is because the 
WIMP-Higgs coupling factor C22 is not suppressed by the mixing angle, but is of the order 
of the standard 4*^ family Majorana neutrino-Higgs couplings, which are already excluded 
by the DM direct detection experiments [40]. However, for Scenario II the WIMP-Higgs 
interaction is suppressed with the mixing angle as C22 — sin^ 9. Thus in Scenario II the 
WIMPs can avoid being ruled out by the spin-independent limits, and actually for this 
scenario the spin-dependent limits give more stringent constraints for the model. 

We start by deriving the constraint following from the SD cross section limits, as these 
were already shown in the previous result section. The best current spin-dependent cross 
section limit comes from the cryogenic dark matter search XENONIO experiment [40]. 
This experiment has given their (spin-dependent) constraints for a standard model 4*^ 
family Majorana neutrino in reference [40] and explicitly plots the expected count rate 
in their detector for this case. Now, the (SD) count rate N oc ao for our WIMP, differs 
from the standard model case only by a simple scaling of the cross section factor ctq: 
c'"o,SM -^ cTQ^Mix = sin^^cTQ^sM, where ctq accounts for the spin-dependent WIMP-nucleus 
cross section at the zero momentum transfer limit. Using this information we can convert 
the XENONIO results for a 4*^ family SM-neutrino to an upper limit on the mixing angle 
as a function of mass: 



The function A/'sm(^^2) was read from the left panel in Fig. 2 of ref. [40] and the function 
Mimit(^^2) was approximated by a linear interpolation between the values of A/'sm(^^2) at 
the high and low mass ends of the SM-exclusion region in the same figure. The white dashed 
line in Figs. 2-4 show the excluded region corresponding to the upper limit (3.1). To make 
this constraint even tighter we project this limit to correspond to the limit following from 
the null results of XENONIOO experiment. This is done by simply scaling the above limit 
for the mixing angle (3.1) with the fourth root of the ratios of the mass-day exposures 
of the two experiments: sin^iim -^ sin^iini(£'xeio/^Xeioo)''^^^- We use the full exposure 
^Xeio — 136kg-days for the XENONIO experiment and exposure £'xeioo — 1471kg-days 
for the XENONIOO, which was reported by the XENONIOO collaboration in connection 
with their April 2011 results in [48]. This tighter limit corresponds to the white solid line 
in Figs. 2- 4. 

As already mentioned above, our WIMP has also spin-independent interactions me- 
diated by the Higgs field, which were not accounted for in the treatment leading to the 
constraint (3.1). The spin-independent limits from XENONlO-100 are given in [39] and 
[48] respectively, in terms of the WIMP-nucleon cross section. In our model this cross 
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Figure 5: Shown are the spin independent WIMP-nucleon cross sections in the case of tuh = 120 
GeV, in Scenario I. pi2 = +1(— 1) in the left (right) plot. Colored patch shows the model predic- 
tions, yielding ^^2 between values 0.19-0.23, and the colormap indicates the value of the invisible 
branching ratio Rr- The black curves line the exclusion regions following from XENONIO [39] 
(dashed), CDMS [41] (dash-dotted), and XENONIOO [48] (solid) experiments. The black dotted 
curve is the full expected XENONIOO sensitivity. For production of the part of the experimental 
curves (XENONIO, CDMS, XENONIOO expected) shown here and in Figs. 6-7 we used the tools 
ofref. [53]. 



section (in the zero momentum transfer limit) is given by 



SI,] 



V^22J Z 114 / 1 



7T 



m 



(3.2) 



H 



where n refers to a nucleon, fi^ is the WIMP-nucleon reduced mass and / is the Higgs 
nucleon coupling factor accounting for the quark scalar currents in the nucleons ^. 

In the figures we have plotted ctq '^ from (3.2) for our model using those 7712, sin0 and 
mi values in C22 which produce the correct DM density Qn^. Thus each point in (ctq '^, 7712) 
plane can be mapped to a point in one of the Figs. 2-4 naturally respecting the tuh and pi2 
cases. Also the exclusion curves following from XENONIO [39] (black dashed), CDMS [41] 
(black dash-dotted), XENONIOO [48] (black solid) results and XENONIOO projected (black 
dotted) are plotted in the figures. The colormap shows the value of the Higgs branching 
ratio Rr to the invisible DM sector. 

4. Conclusions 

It is possible that a nonsequential fourth generation of leptons without quarks exists. We 
have studied the possibility that relatively light Majorana neutrinos may be the dominant 
decay channel of a (composite) light Higgs boson, and that these neutrinos, if stable, may 



'^We use rather conservative value / = 0.5 in our computations. However, the uncertainties in / are 
pretty large, following mostly from the uncertainties in the value of the pion nucleon sigma term (See e.g. 
[49, 50, 51, 52] ). The values of / may vary from ^ 0.3 up to ~ 0.6. 
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Figure 6: Same as figure 5 but for mn = 130 GeV. 
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Figure 7: Same as figure 5 but for mn = 145 GeV. 

constitute the observed dark matter abundance of the universe. However, it seems that 
the area of the parameter space where both these conditions are true is already ruled out 
by direct dark matter detection experiments, the most stringent limits coming from the 
spin-independent dark matter-nucleon cross section. 

We investigated two possible generic realizations of the scalar sector responsible for the 
mass patterns of the fourth generation leptons. In scenario I all masses originate from the 
effective SM-like Higgs field, while in scenario II, in addition to the SM-like Higgs field we 
assume that the mass of the right handed neutrino arises from its interactions with a singlet 
scalar field. In scenario I the area that produces acceptable relic density is practically ruled 
out by the direct dark matter searches. Therefore this model can not produce a realistic 
dark matter candidate. However, to function as an invisible decay channel of the composite 
Higgs boson, the neutrino does not need to be stable on a cosmological timescale. From 
the viewpoint of collider experiments, the effect will be the same as long as the neutrino 



17 




50 60 70 

m_ (GeV) 



50 60 

m. (GeV) 



Figure 8: Same as figure 5 but for Scenario II, with mn = 120 GeV. 
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Figure 9: Same as figure 5 but for Scenario II, witii tuh = 130 GeV. 

is stable on the timescales relevant for the detector. This could be the case even if mixing 
with the SM neutrinos is allowed, provided the mixing angle is small enough. In this case, 
the limits from dark matter searches obviously do not apply, and we conclude that there 
is a large area of parameter space where the neutrinos are the dominant decay channel, as 
shown in figures 2-4. 

In scenario II a large portion of the points that produce the correct relic density lie in 
the area that is not ruled out by the direct dark matter searches. Therefore the neutrino 
in this model, if stable, is a promising dark matter candidate. However, in this part of the 
parameter space the Higgs-neutrino coupling is weak and the invisible decay channel will 
only have a subdominant effect on Higgs searches. The composite Higgs branching ratio 
to the invisible channel in this case is at most of the order of Rr ^0.1. If a light Higgs 
boson is found in the 125 GeV mass range, it will be interesting to study its partial decay 
widths in detail to see if there are any hints of the presence of a subdominant invisible 



18 



0.25 




0.15 



0.05 



50 60 

m_ (GeV) 



50 60 

m. (GeV) 



Figure 10: Same as figure 5 but for Scenario II, with mn = 145 GeV. 

decay channel that would be compatible with this particular dark matter scenario. Also 
this scenario contains areas of the parameter space where the invisible decay channel has a 
large branching ratio, if we assume that the neutrino is stable only on timescales relevant for 
the detector. Thus, if the neutrino-lifetime is suitable, the visible decays of the Higgs boson 
can be highly suppressed in both models (scenarios I and II), as the invisible branching 
ratio is of the order of Rr ^ 0.9 in some parts of the parameter spaces in both models. 

Depending on the value of Rr the composite Higgs could still be found in the standard 
search channels, if the invisible branching ratio is not too large. Then it would just show 
up as having a slightly smaller than expected production rate. On the other hand, if the 
invisible decay channel completely dominates over the SM decays, the Higgs search strategy 
must be altered. In this case the Higgs shows up only as missing energy in the detector, and 
the most promising search channels would be monojets plus missing energy and dileptons 
plus missing energy, as discussed in [54]. One possibility that we have not addressed in this 
work is the case where the neutrino mixes with the SM neutrinos enough to decay rapidly 
inside the detector. If this is the case, and the neutrino pair is the dominant decay channel 
of the Higgs, then this will obviously have significant effects on the Higgs search strategy 
as well. We leave the exploration of these possibilities to future work. 



5. Appendix: Annihilation Cross Section 



Here we give the Majorana WIMP annihilation cross section to SM fermions: 



a 



N2N2^ff 
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27rs /32 ^ 
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2 ^2q2q2 






m 



w 



(5.1) 
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where 9w is the Weinberg angle, (3i = {1 — Amf/s)^/'^ and 

Dx = r-^ ■ (5-2) 

The factor C22 accounts for the different couphng strengths of the hghtest neutral particle 
to the SM-like Higgs as indicated in Table. 1. Finally Nq — 1 for leptons and N^ — 3 for 
quarks and 

gf{s^m2^mf) = {^s{s — mj)/3| + 2im^itnj){v^f + a^) + rnj{s — Qm\){v^^ — aj) 

22 22 

- ( cos^ 9^8^^s - cos^ 9w4.^^s^)a^f , (5.3) 

where Vf = T^f — 2Qf sin^ 9w and a/ = T3/, where T3 is the isospin and Q is the charge of 
the fermion. The cross section in Eq. (5.1) differs from the one given in [12] (and in [36]) 
as the last line in the expression of gf{s,m2,mf) above is absent in [12]. Reason for this 
is that here we use the full Z-boson propagator (in the unitary gauge) when calculating 
the cross section, in contrast to the [12] (and [36]) where only the part of the propagator 
proportional to the ^^^y has been used. This difference has only a minor impact for the DM 
analysis done in this work. The effect of the new terms is strongest for the heavier WIMPs 
(7712 ^ "f^top) foi" which the annihilation channel to top quarks becomes relevant. However 
also in that case, as we are far from the Z-pole, the net effect is not dramatic. 
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